Dynamic spin susceptibility in the t-J model 
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A relaxation-function theory for the dynamic spin susceptibility in the t-J model is presented. 
By a sum-rule-conserving generalized mean-field approximation (GMFA), the two-spin correlation 
functions of arbitrary range, the staggered magnetization, the uniform static susceptibility, and the 
antiferromagnetic correlation length are calculated in a wide region of hole doping and temperaturs. 
A good agreement with available exact diagonalization (ED) data is found. The correlation length 
is in reasonable agreement with neutron-scattering experiments on La2_iSr5Cu04. Going beyond 
the GMFA, the self-energy is calculated in the mode-coupling approximation. The spin dynamics 
at arbitrary frequencies and wave vectors is studied for various temperatures and hole doping. At 
low doping a spin-wave-type behavior is found as in the Heisenberg model, while at higher doping 
a strong damping caused by hole hopping occurs, and a relaxation- type spin dynamics is observed 
in agreement with the ED results. The local spin susceptibility and its ui/T scaling behavior are 
calculated in a reasonable agreement with experimental and ED data. 

PACS numbers: 74.72.-h, 75.10.-b, 75.40.Gb 



I. INTRODUCTION 

It is generally believed that charge-carrier interaction 
with spin fluctuations in the cuprate high-temperature 
Superconductors is the origin of their anomalous nor- 
mal state properties and may be responsible for the 
superconducting transition (see, e.g., Ref. [l|). Inelas- 
tic neutron scattering experiments have revealed quite a 
complicated behavior of the spin-fluctuation spectra in 
cuprates. [3, ^ Therefore, studies of spin fluctuations in 
these materials are essential in elucidating the nature of 
high-temperature superconductivity. Two limiting cases 
can be vi^ell described. In the undoped insulating phase 
the quasi-two-dimensional Heisenberg model for local- 
ized spins gives a reasonable description of the spin- 
fluctuation spectra (see, e.g., Ref. [3|), while in the over- 
doped region the random phase approximation (RPA) for 
weakly correlated itinerant electrons can be applied (see, 
fe.g., Ref. i). 

' However, the region of light and optimal doping (the 
so-called "pseudogap phase"), where localized spins on 
copper sites strongly interact with correlated charge car- 
piers is much more difficult to study. This region should 
be treated within a model of strongly correlated elec- 
trons like the Hubbard model or the t-J model. Q The 
charge-carrier motion in the t-J model is described by the 
Hubbard projection operators, whose commutation rela- 
tions are more complicated than those of Fermi or Bose 
operators. Various approaches have been used to study 
the spin dynamics in the t-J model (for a review see, e.g., 
Refs. [1] and 0). In particular, in the slave boson or 
fermion methods, a local constraint prohibiting a double 
occupancy of any quantum state is difficult to treat rigor- 
ously. An application of a special diagram technique for 
Hubbard operators to the t-J model results in a compli- 
cated analytical expression for the dynamic spin suscep- 



tibility (DSS) . [l^ Studies of finite clusters by numerical 
methods were important in elucidating static and dy- 
namic spin interactions, though they have limited e nerg y 
and momentum resolutions (see, e.g., Refs. [HI, [H, [U ) . 

To overcome this complexity, we apply the projection 
Mori-type technique elaborated for the two-time ther- 
modynamic Green's function (GF). [ij, [l^, [l^ In this 
method an exact representation for the self-energy (or 
polarization operator) can be derived which, when eval- 
uated in the mode-coupling approximation, yields physi- 
cally reasonable results even for strongly interacting sys- 
tems. As our calculations have shown, the decoupling 
of the correlation function of currents, i.e., (dS^/dt), 
in Ref. [l3| is insufficient for obtaining reasonable re- 
sults. Therefore, in the present paper, in studying the 
DSS Y(q,w) = -((S'+|5'Iq))<^ (written in terms of the 
GF |14| ) the mode-coupling approximation in the para- 
magnetic phase is applied to the correlation functions 
of the forces, i.e., {d^S^/dt"^). [l^ A similar approach 
based on the Mori projection technique for the single- 
particle electron GF and spin GF has been used in 
Refs. [18|, il9|, |23j IMI . The magnetic resonance mode ob- 
served in the superconducting state was studied within 
the memory-function approach in Refs. [20l [2^ [23l[2^. 

In this paper we use the spin-rotation-invariant 
relaxation-function theory for the DSS in the t-J model 
derived by us in Ref. [l^ to calculate the static proper- 
ties in the generalized mean-field approximation (GMFA) 
similarly to Ref. (26| and the dynamic spin-fluctuation 
spectra using the mode-coupling approximation for the 
force- force correlation functions. Thereby, we capture 
both the local and itinerant character of charge carri- 
ers in a consistent way. In calculating the static prop- 
erties, in particular the static susceptibility and spin- 
excitation spectrum, we pay particular attention to a 
proper description of antiferromagnetic (AF) short-range 
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order (SRO) and its implications on the spin dynamics. 
For the undoped case described by the Heisenberg model 
our results are similar to those in Refs. [131 and [28l |. 
For a finite doping our theory yields a reasonable agree- 
ment with available exact diagonalization (ED) data and 
neutron scattering experiments. 

The paper is organized as follows. In the next sec- 
tion the t~J model is formulated in terms of the Hub- 
bard operators, and basic formulas for the static spin 
susceptibility and the self-energy within the relaxation- 
function theory [2^ are presented. Numerical results 
for the static properties and spin-fluctuation spectra are 
given in Sec. Ill, where their temperature and doping 
dependences are analyzed. The conclusion is given in 
Sec. IV. Details of the calculations are discussed in the 
Appendix. 



II. RELAXATION-FUNCTION THEORY 
A. Basic formulas 

We start with the t-J model on the square lattice. 



i E J^Axrxr-xrxJ^), (1) 



which is written in terms of the Hubbard operators (HOs) 
X"^ — \i,a){i,(3\ for three possible states at a lat- 
tice site i: for an empty site \i,a) = \i,0) and for 
a singly occupied site \i,a) — \i,a) with spin (7(1/2) 
(cr = ±, a = —a). The HOs obey the multiplica- 
tion rule X°'^Xl^'^ = Xf' and the completeness relation 
Xi^ + So- ^i"^ ~ which preserves rigorously the con- 
straint of no double-occupancy of any lattice site. The 
spin and number operators of the model are given by 
= Xr, Sf = (1/2) <jXr, and = Xr ■ 
The chemical potential ^ is determined from the equation 
for the average electron density n — J2cr(-^i'^) = 1 ^ <5 
where S = (XP°) is the hole concentration. 

In Ref. |25| we have derived the general expression for 
the DSS xici,Lo) = ~{{S+\SZ^))^, 



0J„ 



-cjS(q,tj) 



(2) 



The static spin susceptibility Xq is related to a general- 
ized mean-field spin-excitation spectrum Wq by the equa- 
tion 



Xq = {S^,S_^) = m(q)/t^2 



(3) 



with m(q) — {[iS^ , S_q\) . Here, the Kubo-Mori scalar 
product is defined as (see, e.g., Ref. |16i] ) 



{A{t),B)^ dX{A{t-iX)B), p^l/kBT. (4) 
Jo 



The self-energy is given by [25| 

E(q,c.) = -^M-S+\ - 5lq))£P™P-) , (5) 

where 



iiA\B))^^~t dte''^\A{t),B) (6) 
Jo 

is the Kubo-Mori relaxation function and its "proper" 
part means that it does not contain parts connected by 
a single relaxation function in the GMFA. The spin- 
fluctuation spectrum is given by the imaginary part of 
the DSS (121), 



X"(q,^) 



S"(q, oj) m(q) 



[c.2-,^2_^S'(q,a;)]2 + [u;S"(q,c.)]^ 



(7) 



where S(q,cj + iO+) = S'(q,cj) + iJ:"{q,uj), and 
S'(q,cj) = --S'(q,-cj) and E"(q,w) = S"(q,--w) < 
are the real and imaginary parts of the self-energy, re- 
spectively. 



B. Static properties 

To calculate the static susceptibility and the spin- 
excitation spectrum lu^ in Eq. ([3]), we use the equality 



™(q) = ([z5+,5Zj)-(-^+ ^Iq), 



(8) 



where 



m(q) = -8t{l - 7q)^^i,o - 8J(1 - 7q)Ci,o (9) 

with 7q = (1/2) (cos q^: + cosQy) (we take the lattice 
spacing a to be unity), En^m = Fu = {X^'^ X^), 
Cn,m = C'r {Sq S^), and R = ne^ + mcy . Here, 
we take into account the hopping integral tij and the 
exchange interaction Jij for the nearest neighbors only 
denoted by t and J, respectively. 

To calculate the correlation function (— 5'+,S'Zq) in 
Eq. ([5]), we take the site representation and use the 
decoupling procedure which is equivalent to the mode- 
coupling approximation for the equal-time correlation 
function (see Appendix A). We obtain {—S'^.SZq) ~ 
tjq (5*+ , STq) and, by comparison with Eq. we get 
the spin-excitation spectrum 



= 8t2Ai(l-7q)(l-n-F2^o-2i^i,i) 
+ 4j2(l-7q)(A2^-aiCi,o(47q + l) 
+ a2(2Ci,i+C2,o)). 



(10) 



The decoupling parameters ai,a2,\i, and A2 are ex- 
plained in Appendix A. Thus, the static susceptibility 
can be calculated from Eq. 

The AF correlation length ^ may be calculated by ex- 
panding the static susceptibility in the neighborhood of 



3 



the AF wave vector Q 
^2fc2). We get 

e 



(7i",7r) , XQ+k = Xq/(1 



8J2ai|Ci,o| 



(11) 



The critical behavior of the model ([T|) is reflected by the 
divergence of xq and ^ as T ^ 0, i.e., byti;Q(T = 0) = 0. 
In the phase with AF long-range order (LRO) which, in 
two dimensions, may occur at T = only, the correlation 
function Cr is written as [271, l29l 



(12) 



where Cq = (S^SZ^)- The condensation part C deter- 
mines the staggered magnetization which is defined in 
the spin-rotation-invariant form 



2N 



■ c. 



(13) 



R 



The static susceptibility, the correlation functions, the 
correlation length, and the magnetization are calculated 
in the GMFA for arbitrary temperatures and doping (see 
Sec. III. A). Then, the GMFA results are used for the 
calculation of the self-energy (see Sec. III.B). 



C. Self-energy 

The self-energy ^ can be written in terms of the cor- 
responding time-dependent correlation function as 



S(q,c^) = 



1 



dco'- 



27rrn(q) lu'{uj - uj') 

OO 



(14) 



The self-energy is calculated in the mode-coupling ap- 
proximation for the multisite correlation functions result- 
ing from the operator S^{t) as outlined in Appendix B. 
We consider only the imaginary part of the self-energy 
([T4)) since the real part is given by the dispersion rela- 
tion [13. 

As it turns out by numerical evaluations (see 
Sec. III.B), the largest contributions come from two di- 
agonal terms. For the first term E" we get 



S'i(q,^) 



7r(2 J)4 



-y 



2m(q)LuN(uj) 

\ y qi,q2 



(koi(ko2 



{^qiqaqa + Tqiqsqs ^^qsqiqa l^l'^l )^(^2) (15) 
N{lO -UJl- CJ2)-Bqi {UJI)B^.^{UJ2)B^^{UJ - - UJ2), 

where N{u}) = {e^'^ — 1)^-"^ and q = qi + q2 + • 
The spectral density of the spin-fluctuation spectrum 



Sq(w) = (l/7r)x"(q,w) is given by Eq. The ver- 

tex for the spin-spin scattering reads (cf. Ref. [28j ) 

qiq2q3 = 4(7q3+qi - 7q2)(7 qs 7qi I 

-7qi + 7q3 + 7q2+q3 " 7q2-|-qi • (16) 



The second term E" is given by 



s;'(q,c.) 



i{q)uN{u}) iV2 ^ J_ 



d(jjiduj2 (17) 



{Aqiq2q3 + ^q^q.q J^V (^2 ) "(^ + ^1 - L02)[l - n{ui) 

[(l/4)iVq,(w2) + B^,{uj2)]A^A^i)A^,{lo + - 0-2), 



where n{uj) = (e'^" + 1) 



Here, the single-particle 



spectral function A^{lo) = -(l/7r)Im((X^''|X^"))^ and 
the charge-susceptibility spectral function A^q(w) ~ 
— (l/7r)Im((nq|n_q))(^ are introduced. The vertex for the 
spin-hole scattering reads 



A, 



qiq2 q3 



4(7 



q3+q2 



7qi)7g3 +7q2 -7qi 



-q3 • 



(18) 



The remaining terms in the self-energy are considered in 
Appendix B. In Sec. III.B we calculate the diagonal terms 
for various doping and temperatures. 

We would like to emphasize that in our calculation of 
the self-energy in Eq. (|17p contributions from the charge 
and spin excitations are taken into account explicitly 
by the spectral densities Nq^{L02) and Bq^{L02)- Con- 
trary to this, in Refs. 20] and [22] these terms have 
been approximated by some kind of static or mean-field- 
type expressions. This results in the self-energy of the 
form similar to that given by a conventional particle- 
hole loop diagram used in the weak coupling theory 
like RPA. This form of the self-energy can be readily 
obtained from our expression p7p . if we disregard the 
charge fluctuation contribution and neglect a small spin 
excitation energy lj2 in comparison with the Fermi energy 
in the Fermi function and in the hole spectral function: 
n{u)-\-u)i — L02)Aci3{L0 + uJi—uj2) ^ n{^JJ + ^JJl)Ac^.^{uJ + LOl) . 
Then we can integrate over L02 in Eq. (|17p which gives 
^^^<Lj2N{lo2)B^,{^2) = Cq,, where Cq = {S+SZ^). 
As a result the self-energy takes the form 



S'/(q,^) 



Tl(2tf 1 

m(q)(jj N"^ 



CqalAq^q^q^ -I- Aq^q^q^} 



duJl[n[uj + UJl) - 7l(cJl)]Aq^(a;i)Aq3 [uJ + CJl),(19) 



which is similar to that found in the one-loop par ticle- 
hole approximation used in Refs. [13]- [13] and |22j ] and 
describes the damping due to the decay of spin fluc- 
tuations into electron-hole excitations. The same re- 
sult can be deduced, if in the mode-coupling approxima- 
tion (see Appendix, Eq. (jB2|) the time-dependent spin 
correlation function is approximated by its static value: 
{Xl'^{t)X11^) ^ {X^^XZIJ ■ As our numerical calcu- 
lations have shown, the imaginary part of the self-energy 
(|19p is about twice as large in comparison with that given 



by Eq. (HZD. 
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III. NUMERICAL RESULTS 

To investigate the magnetic properties of the t-J 
model, in particular the spin dynamics at arbitrary tem- 
peratures and hole concentrations, we start from the 
GMFA. By Eq. ^ we get the DSS 



X(")(q,c.) 



2lo„ 



1 



1 



(20) 



with m(q) and Wq given by Eqs. ([9]) and (fTOjl . respec- 
tively, and the correlation function 



m(q) ^^^^ ^ 

2 LOa 2 



(21) 



The electron Green function is calculated in the Hub- 
bard I approximation which yields 



1 - n/2 



with ^ -4(1 - n/2) i 7q. We get 

Fq = (X^°X°-> = (1 - n/2) n{E^ - . 



(22) 



(23) 



The chemical potential jj, is calculated by the number 
condition n = {2/N) J2q F^- 

To go one step beyond the GMFA, we calculate the self- 
energy [Eqs. ([T5l)-([l8l)] by inserting the GMFA results. 
Moreover, for the spectral function Nq{uj) of the dynamic 
charge susceptibility app earing in Eq. P?|) we take the 
GMFA result of Ref. [Sfl]: 

iVW(c^) ^ 8ifi,o(l-7q) 

X il/n^)[Siu;~n^)-S{u; + n^)], (24) 

where il^ ~ 8t^(l — 7q)(l — n/2) in the leading order of 
doping. 



E 0.0 
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FIG. 1; Spin correlation functions vs doping at T = and 
J/t = 0.4 (solid) and the ED data of Ref. [Sli] (symbols). The 
function Ci,o at J/t = 0.2 is plotted by the dashed line. 



limit {S = 0) and their values are taken also at finite 
5. For (5 = we have Fr^o = so that the itiner- 
ant contribution to the spectrum (jTU]) vanishes, and Wq 
agrees with the result of Ref. [131 , where we have to put 
A2 = 1, as can be seen from Eq. (jA9[) . At T = we fix 
ai and a2 by the sum rule Cq.o = 1/2 and, as an input, 
by the value of the nearest-neighbor correlation function 
obtained by exact diagonalization (ED), Cf'^((5 = 0) = 
-0.234 (Ref. [HI). We get ai = 2.285 and = 2.548. 
At finite temperatures we determine ai{T) and a^CT) 
from the sum rule and the ansatz (cf. Refs. [13, H^) 
r(T) = [ai(T) - l]/[a2(T) - 1] = r(0) = 0.830. 

Considering the parameters Ai and A2 at finite S, at 
T = we fix them by the sum rule ([25l) and by the ED 



value Cf^{S = 0.0625; J/i = 0.4) = -0.176 (Ref. '315). 
We get Ai = 0.195 and A2 = 0.515 with A1/A2 0.378. 
At arbitrary temperatures, doping, and ratios of J/t we 
determine Ai and A2 from the sum rule (p5|) and the 

ansatz 



A. Static properties 

Considering first the static magnetic properties in the 
GMFA, we have to solve numerically the coupled sys- 
tem of self-consistency equations for the correlation func- 
tions Cr = (1/^) X]q C'q e*''^ and for the transfer ampli- 
tudes Fr = (1/A^) X]q -^q e"i^. In the long-range ordered 
phase, Eq. and the additional equation ujq — de- 
termining the condensation part C must be taken into 
account. To this end, the parameters ai, a2, Ai , and A2 
have to be determined, where the sum rule 

Co,o = (^o+^o"> = ^ (1 - (25) 

must be fulfilled at arbitrary temperatures and hole dop- 
ing. 

We fix the decoupling parameters as follows. The pa- 
rameters ai and 02 are determined in the Heisenberg 



Xi{T,S; J/t)/X2{T,S; J/t) = 0.378. (26) 

In Fig. [T] our results for the doping dependence of the 
spin correlation functions at T = and J /t = 0.4 are pre- 
sented. They show a good agreement with the ED data 
of Ref. [m. The different signs of C„,™ reflect the AF 
SRO which gradually decreases with increasing doping 
and decreasing ratios J/t. 

Considering the staggered magnetization m{5) at zero 
temperature which is plotted in Fig. ^ we obtain a strong 
suppression of LRO with increasing doping due to the 
spin-hole interaction. In the Heisenberg limit we get 
m(0) = 0.303 which agrees with the value m(0) = 0.3074 
found in quantum Monte Carlo (QMC) simulations [3^ . 
At the critical doping Sc{ J/t) we obtain a transition from 
the LRO phase to a paramagnetic phase with AF SRO. It 
is remarkable that Sc is nearly proportional to J/t. This 
result agrees with that found by the cumulant approach 
of Ref. [33| . where our Sc values are somewhat lower (e. 
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FIG. 2; Staggered magnetization as a function of doping for 
different values of J jt. 



g., in Ref. 'aa^, (5c ~ 0.06 for Jjt = 0.4). Note that the 
proportionahty Sc oc J/t was not found in Ref. [1^ . The 
Sc values obtained are in quahtative agrement with neu- 
tron scattering experiments on La2-5Sr5Cu04 (LSCO) 
which reveal the vanishing of LRO at ~ 0.02. [2|| 




FIG. 3: Uniform static spin susceptibility vs doping at J/t — 
0.3 for various temperatures. The inset shows the position 
5max(T) of the maximum in x vs 5 in comparison with the 
ED data (dots) of Ref. (3|1. . 

In Fig. [3] the uniform static spin susceptibility x — 
(1/2) limg^o Xq at J/t — 0.3 is plotted as a function of 
doping at various temperatures. Within our theory, the 
increase of x(^, T) upon doping is caused by the decrease 
of SRO (cf. Fig.[T]), i.e., of the spin stiffness against orien- 
tation along a homogeneous external magnetic field. At 
large enough doping, x decreases with increasing 5 due 
to the decreasing number of spins. The SRO-induced 
maximum of x a-t Sma.x{T) shifts to lower doping with 
increasing temperature, since SRO effects are less pro- 
nounced at higher T. The doping dependence of x, espe- 
cially the maximum at S^axiT) (see the inset of Fig. E]), 
is in accord with the ED results of Ref. [3J|. Whereas 



the absolute values of x turn out to be lower than the 
ED data, the maximum position is in a remarkably good 
agreement with the ED results. Note that in the spin- 
rotation invariant approach of Ref. [sll for t <^ J , the 
maximum of x as a function of doping was not repro- 
duced. Our results are in qualitative agreement with ex- 
periments on LSCO, where the measured doping depen- 
dence of the magnetic susceptibility exhibits a maximum 
at (5max — 0.25 over the entire accessible temperature 
region, 50 K< T <400 K. [s^ 

Considering the temperature dependence of xC^"; <5) at 
fixed doping, from Fig.[3]it can be seen that there appears 
a maximum at rmax('5) which shifts to lower tempera- 
tures with increa sing doping, in qualitative agreement 
with the ED data [34| . This maximum and the crossover 
to the high-temperature Curie law x c< (1 ~ ^)/T can be 
understood as a SRO effect, in analogy to the explanation 
of the doping dependence. 




T/t 



FIG. 4: Inverse AF correlation length vs T at J/t = 0.4 (solid 
lines) for doping 5 = Q, 0.04, 0.1, from bottom to top, and at 
J/t — 0.2 (dashed line) for 5 = 0.04. The neutron-scattering 
data on La2-i Sr^ Cu04 with 5 = 0.04 are given by symbols 01 • 
The inset exhibits the doping dependence of the correlation 
length at T = (solid line) and T = Q.lt (dashed line) at 
J/t = 0.4. 

Figure 2] shows the inverse correlation length 
^~^{T, (5, J/t). The qualitative behavior of ^ in the zero- 
temperature limit as function of doping and J/t can be 
easily understood by considering the staggered magneti- 
zation at r = depicted in Fig. [51 At a given value 
of J/t and (5 < ^c, in the limit T ^ 0, AF LRO 
emerges which is connected with the closing of the AF 
gap, cjQ 0, and, by Eq. (Ilip . with the divergence 
of ^ . At zero doping, £,^^{T) exhibits the known ex- 
ponential decrease as T ^ 0. [13, [29| ki 5 > Sc, the 
ground state has no AF LRO, i.e., we have ujq > 0, 
and the correlation length saturates at T ^ 0. Equally, 
taking T = 0, the transition from the AF LRO phase 
to a paramagnetic phase with AF SRO a.t S — 5c is 
accompanied by the change £,~^{0,6 < 5c, J/t) = to 
^~^(0, S > 6c, J/t) > . Considering the influence of the 
ratio J/t on the properties of ^, let us compare the curves 
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in Fig. [1 for fixed S = 0.04 and J/t = 0.4 and 0.2. Ac- 
cording to Fig. m for J/t = 0.4 and J/t = 0.2 we have 
(5 < (5c and (5 > Sc, so that C"^(0, 5, J/t = 0.4) = and 
£_^^{0,S,J/t = 0.2) > 0, respectively. The weakening of 
AF correlations (decrease of ^) with decreasing exchange 
interaction is in accord with the results for Ci.o shown 
in Fig. [2 for J/t = 0.4 and 0.2. Note that in Ref. [13] 
a divergence of ^ was found as T — > for arbitrary val- 
ues of (5, in disagreement with experimental facts. This 
deficiency may be due to employing various decoupling 
schemes not used in our theory. 

To compare the temperature dependence of £^~^{T,S) 
with neutron-scattering experiments on LSCO at T < 
600 K, [2] we take a = 3.79 A and J = 130 meV and con- 
sider the doping S — 0.04. As can be seen in Fig.|4l we ob- 
tain a reasonable agreement with experiments. Concern- 
ing the doping dependence of ^(J, T) depicted in the inset 
of Fig. m it can be described approximately by the pro- 
portionality ^((5, T) cx l/VS (dashed line) which agrees 
with the experimental findings. 0] 



B. Spin dynamics 

In this section we present results for the spin- 
fiuctuation spectra provided by the imaginary part of the 
DSS x"(q, w), Eq. ([7]), where we neglect the real part of 
the self-energy I]'(q, w) (cf. Ref. [2^). The damping of 
spin fluctuations r(q, uj) is determined by the imaginary 
part of the self-energy, r(q, w) = — (l/2)I]"(q, a;) , con- 
sidered in Sec. II. C. Here we mainly consider the damping 
aX UJ — Uq, Fq = r(q, Ljq) . It turns out that the major 
contributions to the damping are given by the diagonal 
terms Ej(q,u;) , Eq. ([15]), and S"(q,a;) , Eq. HT]), while 



the interference terms, such as S'Jj j((q, w) , Eq. (jB4[) . 
appear to be much smaller and may be neglected. That 
is, the damping Fq is the sum of the spin-spin scattering 
contribution Fj_q = — (l/2)Ej(q, Wq) and the spin-hole 
scattering contribution Ft^q = — (l/2)S^(q, Wq) 




(0, 0) 



0) 



{n. k) 



(0, 0) 



FIG. 5: Spectrum ujq (solid line) and damping Fq (dashed 
line) in the Heisenberg limit, 5 = 0, at T = 0.35J . 
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FIG. 6: Dynamic structure factor ^(q,^;) in the Heisenberg 
limit, 5 = 0, at T = 0.38J for wave vectors: qi — 7r(l/2, 1/2) 
(solid line) and q2 = 7r(5/8, 5/8) (dashed line) in comparison 
with the QMC results of Ref. [39( ] for qi (dashed-dotted line) 
and for q2 (dotted line). 



rj,q + Ft,q . Note that in Refs. ^ and ^ the partition of the damping to the spin- wave excitation energy tends 



of the damping into a spin-exchange contribution and a 
fermionic contribution was suggested from the ED data. 
The numerical calculations of S"(q, cj) are performed 
for the exchange interaction J = 0.3i, the value which is 
usually used in numerical simulations. This affords us to 
compare our analytical results with finite cluster calcula- 
tions and to check the reliability of our approximations. 

Let us first consider the Heisenberg limit S — 0. Fig- 
ure [5] shows the spin-excitation spectrum Wq, Eq. ()10|) . 
and the damping Fq = Fj.q. The results are similar to 
those obtained in Ref. I 



In the spin- wave region, at 
^ 1, we get well-defined quasiparticles with Fq ^ cjq. 
For example, for q = tt (1/2, 1/2) and T = 0.35J we 
have g > ~ 0.1 (see Fig.© and Fq/cjq ~ 0.1. Well- 
defined spin excitations for the two-dimensional Heisen- 
berg model have been found by several authors (for a 
review see Ref. Q). In particular, as shown in Ref. [s^ . 
if T ^ and q ^ with the restriction q£, ^ 1, the ratio 



to zero: Fj^q/wq —i- 0. 

To compare our results for the damping with the QMC 
data of Ref. we have considered the linewidth Aq 
of the relaxation function F{q,uj) = A\0uj'Xa]~^'x" (q, m) 
at T = 0.35 J, where Aq ~ 2Fq (Ref. [2^), and have 
found a good agreement. For a further comparison with 
the QMC data we calculate the dynamic structure fac- 
tor S{c[,uj) = [1 — exp(— /3a;]~-^ x"(qi ^)- In Fig. [5] our 
results at T = 0.38J and for two wave vectors in the 
spin-wave region are plotted. The peaks in 5(q, cj), oc- 
curring nearly at Wq (cf. Fig. [5]), reveal well-defined spin 
excitations. Comparing the peak heights with the QMC 
values, we get a better agreement than it was found in 
Ref. [28| . This may be ascribed to Eq. (fT5| which corrects 
the result of Ref. [2^| by the appearance of the additional 
term Fqjq^qgFqjq^qj that cannot be written as a square 
of the vertex Fq^q^q,. 

For non-zero doping the spin-hole scattering contri- 
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FIG. 7: Spectrum ujq (solid line), and damping Fj.q (dotted 
line) and Ft^q (dashed line) at T = 0.15t and 5 = 0.1. 
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FIG. 8: Doping dependence of damping Fq at T = O.lt (a) 
and T = 0.15t (b). 



bution E"(q, w), Eq. ([T7)) . increases rapidly with dop- 
ing and temperature and already at moderate hole con- 
centration far exceeds the spin-spin scattering contri- 
bution Ej(q, a;), Eq. (fT5|) . as demonstrated in Figs. [7] 
and [51 Depending on q, doping, and temperature, the 
spin excitations may have a different character and dy- 
namics. In particular, for the spin-spin scattering con- 
tribution Ej(q, w) we observe, in the long- wavelength 
limit, limq^oTj^q — 0, as in the case of the Heisen- 
berg model shown in Fig. [5l Contrary to this behav- 
ior, the damping F^ q, induced by the spin-hole scat- 
tering, is finite in this limit, both taking Eq. (fTTI) and 
Eq. ([TO)) . The different behavior of Fj.q and Ft.q may 



be explained by the different q-dependence of the spec- 
tral functions entering Eqs. (|15p and p?)) . Whereas for 
spin excitations the spectral function is proportional to 
m(q)/a;q ^ g for g ^ (see Eq. (I^Hl) ). for electrons it is 
finite in this limit (sec Eq. ((22)) ). Therefore, in the limit 
of q = qi + q2 + qs = , for the spin-spin scattering 
the product m(qi)m(q2)TO(q3)/wqjCjq2u;q3 gives a van- 
ishingly small contribution to the integrals over qi , q2 
in Ej(q, w) , Eq. (fTS)) . while in the spin-hole self-energy 
(fT7)) there is no such a small factor. The physical mean- 
ing of the finite damping Ft q at g = can be explained 
similarly to that of the finite electrical conductivity in 
the low- frequency limit, which is in fact a response func- 
tion determining the damping of charge fluctuations at 
q = 0. As is well known, the relaxation rate for the con- 
ductivity at zero frequency, i.e., the inverse resistivity, is 
finite, if one takes into account momentum relaxation of 
electron-hole pairs on phonons. 
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FIG. 9: Damping Fq = -(1/2) E"(Q,a; = u;^) (a) and low 
energy damping F(Q) = -(1/2) S"(Q,c^ = 0) (b) as func- 
tions of temperature and doping in comparison with the ED 
dataofRef. ^ (symbols). 

To discuss the temperature and doping dependence of 
the damping (see Fig. [5]) in more detail, we choose q = Q 
and consider the damping Fq = — (1/2) S"(Q, cj^) as 
function of T and i5 that is plotted in Fig. [9] (a). In the 
zero-temperature limit and for S < d^, where there is 
AF LRO (see Fig. (2), it can be shown analytically that 
Fq(T = 0,5 < 6c) = 0. That is, in the LRO phase we 
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get well-defined spin waves. The vanishing of Fq may 
be explained as follows. Spin excitations at T = can 
decay into particle-hole excitations with a positive en- 
ergy Lo only to satisfy the energy conservation law. Here, 
Lo = uj^ — for 5 < 6c so that we have no damping. On 
the other hand, at T = and 5 > 6c, there is no LRO 
and > which resuhs in FQ(r = 0,6 > 6c) > (see 
Fig. [S] (a) for 5 = 0.1 and 0.15). With increasing tem- 
perature and doping the damping increases as expected. 
To compare our results with the data of Ref. [21|, in 
Fig. [5] (b) the temperature dependence of the low-energy 
damping F(Q) = —(1/2) 'E"(Q, uj = 0) for various dop- 
ing is shown. In Ref. [2l| the function 7(Q) = 2r(Q) 
was extracted from the finite-T Lanczos data for the 
spectral function x"(g,a;) using a simplified ansatz for 
the spin-excitation spectrum, Wq cx ((q — Q)^ -I- k^) with 
K = taken as a temperature- independent parameter, 
whereas our theory allows a direct microscopic calcula- 
tion of F(Q). In the high-temperature region a remark- 
ably good agreement is found. In the low-temperature 
region, which is not accessible by the finite-T Lanczos 
method, the ED results were extrapolated to T = with 
a finite value ofF(Q;T = 0,i5 > 0). This is in contrast 
to our resuh F(Q;T = 0, (5) = (cf. Fig. [9] (b)) which 
agrees with rQ(r = 0,6 < 6c) ^ (cf. Fig. E (a)), be- 
cause both quantities are calculated for to — 0, and which 
may be understood as explained above. 

As illustrated in Fig. [5] and Fig. [51 at low enough 
doping and temperature, i.e., at small enough Tt,q, we 
may observe well-defined high-energy spin-wave-like ex- 
citations with q,k^ 1/^ (fc = |q — Q|) and Fq ^ Wq 
propagating in AF SRO. Considering, for example, spin 
excitations with q ~ (7r,0) at J = 0.1 and T = 0.15^, 
we have qS, — 5.8 with ^ — 1.85 taken from Fig. HI 
Wq = 0.66* , and Fq = 0.27* (see Figs.[7]and|H](b)). That 
is, in this case we have strongly damped spin waves. 

To discuss quantitatively the spectral function 
x"(q, uj) shown in Figs. [TOl and [TTl in particular the posi- 
tion of its maximum at ujm , we first simplify Eq. ([7|) . By 
numerical evaluations, we have found that the imaginary 
part of the self-energy only weakly depends on frequency, 
which qualitatively agrees with the results of Ref. [2l| . 



-'/'7q < ^q/^q ^ 1 and using Eq. ([3]) we get 



Therefore, we put — E"(q, w) = -q^ 
Eq. ([7|) we get the resonance form 

,2\2 



2rq. Then, by 



X"(q,w) = m(q) — 



(27) 



which has a maximum at uj^ given by 



0.^ = ^{2u.l -4+ [12c.^ + {2u;l - vDT'V^', (28) 

where lim^^o^i^m = ^q- 

Let us consider the region of low-frequency over- 
damped spin-fluctuation modes playing an important 
role in the normal phase of the cuprate superconductors, 
i.e., uj < LOq <^ rjq . Expanding Eq. (j27p with respect to 



X"(q,^) = XqTq 



— ^; fq^^, (29) 



where Fq is the spin-fluctuation excitation energy. Con- 
trary to Ref. [2l|, where a similar expression was de- 
rived, we do not use an ansatz for the spin-excitation 
spectrum ujq (see above) , but calculate it microscopically 
by the GMFA [see Eq. Unj)]. The Lorentzian ^ has 



a maximum at w,^ = Fq, which may be also obtained 
from the expansion of Eq. ((28l) with respect to (wq/yyq)^. 
The overdamped form corresponds to the susceptibility 
x(q, w) = Xqrq(rq — zw)"^ and, as a phenomenologi- 
cal ansatz, has been frequently invoked in the study of 
cuprates, e.g., in the calculation of the normal-state spin- 
fluctuation conductivity. |40|] 




0.5 „/t 1.0 



FIG. 10: Spectral function x"{^,^) for q = Q = (tt,?!") at 
T = O.lt (a) and T = 0.15t (b) for S = 0.03 (solid line), 
S — 0.06 (dashed line), and 5 = 0.1 (dotted line), and for 
q = (vr, 0) at r = O.lt (c) [note the change in the energy 
scale] . 

To exemplify the spin dynamics in different regions, 
we first consider spin waves with q = (tt, 0) at 6 — 0.1 
and T = 0.15t (see above). We get (see Fig. [TTl) 
uJm = 0.59t =i = 0.60* and = 0.81* . That is, 
we find those excitations to have a resonance character. 
On the other hand, the ED data of Ref. [2l[ yield evi- 
dence for an overdoped spin dynamics. This difference 
may be due to a slight underestimation of the damping 
in our theory, which can be also seen in the more pro- 
nounced peaks of the dynamic structure factor at 6 = 
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FIG. 11: Spectral function x"{<Ij^) for various wave vectors 
at T = 0.15t and 5 = 0.1 in comparison with ED data (filled 
symbols, Ref. [U). 

(see Fig. [S]) as compared with the QMC data. As seen in 
Fig. [TO] (c), with increasing S the maximum in x"(q, w) 
with q = (tt, 0) is shifted to lower frequencies, in quah- 
tative agreement with the theory of Ref. ^is'l . 

Next we consider the spectral function at q = Q. At 
very low doping, e.g., S ~ 0.03, and low enough tem- 
perature the damping Fq is very small (see Fig. |S] and 
Fig. [9] (a)), where Fq ^ cjq. In this case we observe 
underdamped spin modes characterized by sharp reso- 
nance peaks in x"(Q,^), as seen in Fig' [TO] (a), (b). 
With increasing doping those modes evolve into over- 
damped (relaxation-type) spin-fluctuation modes (AF 
paramagnons) described by the broad spectrum ([^5)1 . 
For example, considering the AF mode at 8 = 0.1 and 
T = O.lbt (see Figs. [3 [S] (b), [TO] (b), and [11]), we have 
ujQ = OAt, Fq = 0.24t, and = 0.32t ~ cof^ = 0.33t. 
That is, the spectrum of this mode may be well described 
by the overdamped form (|^^ . As seen in Fig. [TT] our 
spin-fluctuation spectrum is in a remarkably good agree- 
ment with the ED data of Ref. [2l[ . Let us consider the 
shift of the maximum in x"(Q,a;) at ujm with increas- 
ing doping at fixed temperature. As can be seen from 
Fig. [TO] (a), (b), at low (high) temperatures, uJm slightly 
increases (decreases) with doping, which results from the 
doping dependence of and Fq. The increase of tOm 
with S at low T is in qualitative agreement with the find- 
ings of Ref. [l^ (T = 0.02t) and with experiments. 

In Fig. [TO] the dynamic structure factor 5(q, w), re- 
sulting from the spectral function shown in Fig. [TTJ is 
plotted. At cj = 0, by Eq. (O, we have 5(q, 0) = 
rxqS"(q, 0)/Wq, and for overdamped modes [Eq. (|29)) ] 

we get 5'(q, 0) = Txq/Fq- The shape of S'(q, uj) for para- 
magnons is in a marked contrast to that for spin waves 
(compare also with Fig. [6]). 
Finally, we present results for the local susceptibility 

XLM = ^Ex"(q>^)> (30) 

q 



4 




FIG. 12: Dynamic structure factor S{q,iL!) a,t T — 0.15t and 
S = 0.1 for various wave vectors. 
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FIG. 13: Local spin susceptibility x'Li'-^) a.t T — and 
small doping. The intensity scaled experimental data on 
Lai.96Sro.o4Cu04 at T = 10 K are shown by dots where 
t = 420 meV is taken. 



by using the data for x"(q,t^), Eq. dZj). In Fig. [TO] the 
local susceptibility at T = and small doping, 6 = 0.04 — 
0.051, is shown. In the neighborhood of the AF phase 
transition at T = and S = Sc === 0.037 (see Fig. [2 
luq = 0), the spin excitations with q ~ Q are weakly 
damped (see Fig. [S] (a)). Therefore, at sufficiently low 
6, the local susceptibility, which is just the density of 
states (DOS) for spin-fiuctuations, reveals a resonance 
maximum at a frequency being close to loq with a high 
DOS (see Fig. [TO] at (5 = 0.05). Because loq and Fq 
decrease with decreasing 5, the maximum shifts to lower 
frequencies and becomes very sharp (in Fig. [TO] at S ^ 
0.04, only the upturn with decreasing frequency is seen). 
On the other hand, with increasing doping the damping 
becomes large enough to wash away the maximum (cf. 
Fig. [TO]) . Note that we obtain, in addition to the low- 
energy maximum, a broad maximum in x'Li^) 
maximum energy of spin excitations, lj 2J = 0.6t (cf. 
Fig. [5]) . This feature was not found in Ref. [2l[ , since a 
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simplified spin-excitation spectrum Wq was used. 

1.2 I , , , , 
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FIG. 14: Scaling function f{uj/T) for various temperatures at 
doping S — 0.04. The solid line is the scaling function found 
in the neutron-scattering experiments on Lai.96Sro.o4Cu04 
(Ref. [HI) and given by Eq. pT)!. 

The pronounced upturn behavior is observed in 
neutron-scattering experiments on lightly doped cuprate 
compounds at low energies and temperatures (see, e.g., 
Refs. [Ill and We get a reasonable agreement 

with experimental data for Lai.96Sro.o4Cu04 (w < 
50 meV) [4l|, if we take the energy scale t = 420 meV 
(see Fig. [T31) , which is the standard value of t ~ 400 meV 
in the t-J model for cuprates. A much better agree- 
ment with the data in Ref. |4l| can be obtained at the 
finite temperature T = O.Olt but for the energy scale 
t = 1.1 eV. A qualitatively similar behavior has been 
found in Ref. |21j within a semi-phenomenological the- 
ory, where the agreement with experiment was achieved 
by the choice i ~ 0.1 eV. 

Figure [U shows the scaling function f{uj/T) = 
xl{uJ,T)/xliuJ,T = 0). The scaling behavior is in 
a remarkable agreement with the data of the neutron- 
scattering experiments on Lai.96Sro.o4Cu04 (Ref. filj ) 
which is shown by the solid line and described by the 
function 



/ (tp;^ = ~ arctan 



ai 



T 



02 



T 



(31) 



with ai — 0.43 and a2 = 10.5. A similar scaling was ob- 
served in the underdoped YBa2Cu30fi,35 with the param- 
eters oi = 0.9 and 02 — 2.8 (Ref. (i^l). Our results may 
be well approximated by the scaling function (pij) with 
ai = 3, but without the (w/T)'^ term, 02 — 0, that gives 
a nonlinear behavior at low values of {i^/T) <C 1. The 
weak non-monotonous behavior of our scaling function at 
T = O.OH results from the appearance of a flat maximum 
in Xl('^:2^) (see above). Note that in Ref. [llj the same 
scaling function (pij) was found with ai = 1.2, 02 = , 
that results in the saturation /(w/T) ^ 1 at higher val- 
ues of {u> /T) > 2 and in strong deviations from the exper- 
iments on LSCO reported in Ref. [4l[ , but yields a good 
fit to experiments on Zn-substituted YBa2Cu306.6 p3| . 



This variation of the scaling function may be explained 
by different values of doping and of the corresponding pa- 
rameters determining the scaling behavior (e.g., the AF 
correlation length ^, Ref. '2l'|). 

To sum up, our studies of the DSS show a crossover 
from well-defined spin-wave-like excitations at low dop- 
ing and temperatures to relaxation-type spin-fiuctuation 
excitations with increasing hole doping, which is in agree- 
ment with inelastic neutron-scattering experiments and 
numerical simulations for finite clusters. Moreover, we 
observe a remarkable agreement of the scaling func- 
tion with the data of neutron-scattering experiments on 

LSCO. m 



IV. CONCLUSION 

The relaxation-function theory for the DSS in the t-J 
model in terms of Hubbard operators is formulated. By 
using a spin-rotation-invariant theory for the DSS derived 
by us in Ref. [2^, we calculate the static properties in 
the GMFA similarly to Ref. [2^ and the spin-fluctuation 
spectra. The mode-coupling approximation for the force- 
force correlation functions, which take into account both 
the exchange and kinetic contributions, was used. For 
the undoped case described by the Heisenberg model our 
results are similar to those in Ref. [28=| and for finite dop- 
ing they show a reasonable agreement with available ED 
data and neutron scattering experiments. 

Contrary to the previous studies based on the memory- 
function method in Refs. [l^ - [2l[, we have taken into 
account all contributions to the spin-excitation spectrum 
Wq in the GMFA and to the self-energy S(q, oj) and thor- 
oughly analyzed their temperature and doping depen- 
dence. In particular, we have found that the contribu- 
tion from the hole- hopping term cx in the spectrum oiq, 
Eq. pH)) . is large even in the underdoped region, 5 < 0.1, 
and results in a rapidly increasing with doping gap ljq at 
the AF wave vector Q. This increase is much larger than 
in the calculations in Ref. [3 1 where the contribution 
has not been considered. We have also shown that the 
largest contribution to the self-energy S(q, oj) comes from 
the hole-hopping term I]"(q, cj) oc t'*, Eq. pT|). at finite 
doping (see Fig. [7]). This is in accord with Ref. [13], but 
contrary to the approximation in Refs. [ll| , [l^ , and [13] , 
where only the mixed contribution E'Jj jj(q, w) cx J'^t^ 
has been taken into account. The latter approximation 
should strongly underestimate the damping of spin exci- 
tations at finite doping. In our calculations of E"(q, w) 
we have also considered explicitly a contribution from 
spin excitations in the self-energy (see Eq. pT)) ). while in 
Refs. [20*] and [23] this contribution has been considered 
in some kind of static or mean-field- type approximations. 

A comparison of the DSS derived within the memory- 
function approach, Eq. ([7]), with the RPA form x(q, w) — 
Xo(q,w)/Jl - gqXo(q,t^)] (see, e.g., Ref. [5]) has shown 
that the RPA expression provides reasonable results at 
large doping, while in the underdoped region the RPA 
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formula fails to describe spin-wave-like excitations. [2j] 
Whereas the damping of spin excitations in Eq. ([7]) at 
low doping is quite small, e.g. Fq 0.2t at (5 = 0.1 (see 
Fig. [7]), within the RPA it is much larger, T ^ t. This re- 
sults in ovcrdamped spin dynamics described by Eq. (|29p 
even in the underdoped region. Thus, we conclude that 
the relaxation-function approach is a reliable theory for 
studying the spin dynamics in a broad region of doping 
and temperatures. 

In this paper we have not performed a fully 
self-consistent calculation of the electronic and spin- 
fluctuation spectra by evaluating the spin correlation 
functions ((2T|) in the GMFA and the electron correla- 
tion functions (I23p in the Hubbard I approximation. As 
was shown in Ref. [31, static AF spin correlations and 
self-energy effects result in a strong renormalization of 
the electronic spectra and should be taken into account 
in a consistent theory. This generalization will be con- 
sidered in a subsequent publication. The theory will be 
also formulated for the superconducting state by intro- 
ducing matrix electronic GF with normal and anomalous 
components as given in Ref. [i^ . 
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APPENDIX A: DECOUPLING PROCEDURE 

To calculate the correlation function (— Sj^, S*;^) in 
Eq. ([H]) we consider the equation 

- S+ = [[S+, {Ht + Hj)] , {Ht + Hj)] EE J2 Ft, (Al) 

a 

where Ht and Hj are the hopping and the exchange parts 
of the Hamiltonian ([T]) and a = tt, tJ, Jt, J J. Here we 
have 

- E {^^" + - ^ } , (A2) 

= It. '^^4'^^" i^^*^" + "'^"1 - ^ .?)}, (A3) 

Ff' = [[S+,Hj],Htl = [[S+,Ht],HjUAA) 

where 

p qz qz q+ , q+ qz qz 

^ijn — ' n i j 

— — S^^S^ , (A6) 



Hy „ — S^Sj + S^S^ Sj 

— S^Sj'Sj^ — S'„ S^Sj' . (A7) 

Explicit expressions for Fl'^Ff* are given in Ref. [25l |. 

To evaluate the corresponding multiparticle correla- 
tion functions in (— 5q,5Zq) we perform the follow- 
ing decouplin g p rocedure similar to that proposed in 
Refs. [26|, l28l. l29l Isst preserving the local correlations. 
The correlation functions from H^^^ are decoupled as 

(xf"x+-xr , sn = Al (xf o^r) {s+, s-), im 

where, for n = i, {X[°X^'') = (Zf ) = n/2 and 
the second term of with n i is neglected (cf. 

Ref. [1^). Decoupling the correlation functions from 
Tlijn=i and Pij^n=i we introduce the parameter A2: 

(H,,„ s^) = -{[X++ + X--] s+, SD 

^-X2{[X++ + Xr-]){S+,S^); 

{p^,^,sr) = -(1/2) i[x++ + xr-]s+, sn 

= -A2 (1/2) {[X++ + Xr-]) {S+, sn, (A9) 

where we used the equations: S^S^ = X^^ , SfS^ + 
S+Sf = 0, and {Sff = (1/4)(X++ + Xr") . Here 
we take A2 7^ Ai , in contrast to the approach of 
Refs. m, li^, A2 = Al . In the Heisenberg limit 5 — we 
have Xf^ + X^ = 1 so that A2 = 1 . The parameters 
Al, A2 describe the renormalization of the vertex for spin 
scattering on charge fluctuations. 

Considering the correlation functions from Hy„^i and 
Pijn^i, where {i, j, n} forms a nearest-neighbor sequence, 
we apply the decoupling scheme used in Refs. [26l. [28l. : 

{S+S+S-,Sn = a,{S+S-)iSl,Sn 

+ a2{S+S-)iS+,S^). (AlO) 

Here, the parameters ai and a2 are attached to nearest- 
neighbor and further-distant correlation functions, re- 
spectively, and describe the renormalization of the vertex 
for spin-spin scattering. The determination of all param- 
eters is considered in Sec. HI. A. 

Calculating the spin-excitation spectrum and the static 
susceptibility ([3]) with Eq. ([8]) we may take into ac- 
count only the diagonal contributions (jA2p and (|A3p 
and omit the mixed contributions corresponding to Ff* 
and Fl^^ in Eq. (|A4[) according to the following rea- 
soning. The mixed contribution of the type (F/*, Sf) 
in the GMFA is proportional to the difference of the 
correlation functions of the form {X^^ — X^^) or 
(A:,°"X°" - , which vanishes in the param- 

agnetic phase. In the same approximation the mixed 
contribution of the type (-F/"^, 5;") turns out to be pro- 
portional to higher-order correlation functions of the type 
(X+"X°", 5";") and may be neglected. 
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APPENDIX B: MODE-COUPLING 
APPROXIMATION 



Using Eq. (lAll) the two-time correlation func- 
tion in Eq. (|14p yields 16 terms for the self-energy 
S(q,w) = J2a,f:s^a,i3{(i,uj), where, e.g., Ejxj,7(q,w) = 

((-^q In the site representation of 5+ 
in Eq. given by Eqs. (|XT|) - (|T7)) we take into ac- 
count products of three spin operators on different sites 
only. This is clear in the Heisenberg limit [2^, where 
terms with coinciding sites reduce to single operators and 
the proper part of the correlation function in Hjjjj is 
considered. For finite doping, as revealed by numerical 
evaluations, the exclusion of terms with coinciding sites 
yields a better agreement with exact data than the inclu- 
sion of those terms. We calculate the two-time correla- 
tion functions in the mode-coupling approximation (see, 
e.g., Ref. |28|), i.e., we approximate them by a product 
of three single-particle two-time correlation functions as 
follows: 



X (<5ki,-k'i^k2,-k^ + '5ki,-k^'5k2,-k'i) (^ka,- 



(Bl) 



(B2) 



In Ref. [25| we have shown that in the Born approxi- 
mation (i.e., in the second order of the effective vertices 
^^ J2,tj) only six contributions to the self-energy may 
be retained, so that 

E(q, = T,jjjj{c[,uj) + Stt,(t(q, w) -1- Stj,tj(q, 

+ I]jt,jt(q,w) + 2Stj,jt(q,c^). (B3) 

The imaginary parts of the diagonal terms Hjj^jj = Sj 
and Yjtt.tt = are given by Eqs. (fT5|) and (fT7| . For one 
of the interference terms we obtain 



qi.q2 



(koi(hj2N {ijji)n{uj2)n{uj — uji — loi) 



(B4) 



with Fqjq^qg glvcu by Eq. (|16p . where the contributions 
linear in 7q reflect the exclusion of terms in S'^ with 
coinciding sites. 
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